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ABSTRACT 

New high energy emission features have been recently discovered by the Chcrcnkov 
telescopes from active galaxies e.g., a few minutes variability time scale of TeV emission 
from Mrk 501 and PKS 2155-304, sub- TeV 7-ray emission from GeV peaked blazar 
3C 279, and TeV emission from two nearby active galaxies, M87 and Cen A, which 
jets are inclined at a relatively large angle to the line of sight. These results have put 
a new light on the high energy processes occurring in central parts of active galaxies 
stimulating more detailed studies of 7-ray emission models. Here we report the results 
of a detailed analysis concerning the most general version of the model for the 7-ray 
production by leptons injected in the jet which interact with the thermal radiation 
from an accretion disk (the so called external inverse Compton model). We investigate 
the 7-ray spectra produced in an anisotropic Inverse Compton (IC) e ± pair cascade in 
the whole volume above the accretion disk. The cascade 7-ray spectra are obtained for 
different locations of the observer in respect to the direction of the jet. We also study 
the time evolution of this 7-ray emission caused by the propagation of the relativistic 
leptons along the jet and the delays resulting from different places of the origin of 
7-rays above the accretion disk. We discuss the main features of such a cascade model 
assuming constant injection rate of electrons along the jet. We are investigating two 
models for maximum energies of injected electrons: with a constant value independent 
on the distance along the jet or limited by the synchrotron energy losses considered 
locally in the jet. The model is discussed in the context of blazars observed at small 
and large inclination angles taking as an example the parameters of the two famous 
sources Cen A and 3C 279. 

Key words: galaxies: active: individual: Cen A, 3C 279 — radiation mechanisms: 
non-thermal — gamma-rays: theory 



1 INTRODUCTION 

Recent observations of active galaxies in the TeV 7-rays con- 
firmed earlier reports on an extremely variable nature of this 
emission. It can c hange on various tim e scales starting from a 
few minutes (e.g . iGaidos et alj j 19961 ): I Albert et alfll2007t): 
Aharonian et al, ( 20071 ) ) up to months/years (|Hsu et all 
20091 ). This variable emission can correspond to the time 
scales of, either the periods of enhanced energy generation 
in the disk-jet system, or acceleration time of particles, or 
the passing time of relativistic particles through the medium 
which creates the target for efficient 7-ray production. More- 
over, the time scales can be modified due to relativistic ef- 
fects caused by the movement of the emission region to- 
wards the observer and, for very distant sources, cosmolog- 
ical propagation. The precise location of the emission site 
within the active nuclei remains unclear. Theoretical argu- 
ments suggest that the emission place is rather distant from 



the black hole. 7-ray radiation produced closer cannot es- 
cape to the observer due to a very strong, soft radiation 
field originating from the inner part of the accretion disk. 
On the other hand, this strong radiation field can create 
a target on which 7-rays could be produced. Recent ob- 
servations of the nearby active galaxy, M87, whose jet is 
observed at a relatively large angle, suggest that the TeV 
7-ray emission region can likely appear wi thin a few tens of 
Schwarzschild radius from the black ho le ( Aharonian et all 
120061 : lAcciari et all 120081 ; lAlbert et all l2008al ). It is likely 
that the emission region is close to the accretion disk inside 
the inner jet as suggested by the recently discovered connec- 
tion between the TeV 7 -ray flare and the su bsequent appear- 
ance of the radio blobs jAcciari et al ] |2009l ). Also Cen A, an- 
other nearby active galaxy inclined at larg e angle, has been 
recently detected as a TeV 7-ray source (|Aharonian et all 
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l2009h . This emission is as well i nterpreted as originated ve ry 
close to the central black hole (jRieger fc Aharonianll2009l ). 

Recently the LAT Telescope on board of the Fermi 
observatory is operating with an order of magnitude bet- 
ter sensitivity than the EGRET telescope. So far, after 11 
month s of operations it observed 709 AGNs in the GeV 
range (|Abdo et alJl201Cft . Among them 300 are BL Lac ob- 
jects and 296 are flat-spectrum radio quasars. On top of this, 
8 objects have been associated as misaligned blazars. There- 
fore, investigation of the anisotropic high energy radiation 
processes in the active galaxies is very timely. 

It is widely believed that 7-rays are produced in active 
galaxies by particles (leptons or hadrons) accelerated in the 
jet launched perpendicularly to the plane of the accretion 
disk. The acceleration mechanism is not well known. Usually 
two idealized cases of angular distribution of particles in 
the jet are discussed, either isotropic in the jet frame or in 
the recti-linear acceleration regions generally directed along 
the jet. The first case corresponds to the acceleration of 
particles in a shock wave. Due to a mild relativistic motion 
of the emission region this scenario is sometimes referred 
to as a slow jet model. The latter case may correspond to 
the acceleration in magnetic field reconnection regions. This 
scenario is sometimes called as a fast jet model because of 
an ordered, super-relativistic motion of particles. 

7-rays can be produced in different radiation mecha- 
nisms. A few basic scenarios have been proposed soon after 
discovery of the GeV 7-ray emission from blazars by the 
EGRET telescope on the CGRO. Among the famous lep- 
tonic models are: syn chrotron self-Compton (SSC) model 
l|Maraschi et al.lll992l ). external Compton (EC) mod el with 
scattering of the disk ra diation (Dermer et al.lll992l ). or ra- 
diation around the disk (|Sikora et al.l 1994 ). Also hadronic 
models have been discussed, e.g. proton in itiated cascade 
(PIC) model (|Mannheim fc Biermannlll992l ). Starting from 
that discovery period, several modifications of these first 
more general scenarios have been discussed in detail with 
the applications to specific objects. A more realistic three 
dimensional cascade model occurring in the strong radiation 
of the accretion disk are starting to be considered only re- 



centl y (|Sitarek fc Bednarekl2010l ; lRoustazadeh fc Boettcherl 
120101 ). In this paper we are specially interested in the for- 
mation of the 7-ray spectrum in terms of the external 
Compton accretion disk mod el (see for details of this model 
iDermer fc Schlickeiserl (| 19931 )). 

It has been shown that for typical parameters of blazars 
detected in GeV-TeV 7-ray energies, the optical depths for 
7-rays injected relatively close to the disk at an arbitrary 
place can be clearly abov e unity, see e.g. recent calcu lations 
for 3C 279 and 3C 273 (ISitarek fc Bednarekl feoOS) or for 
Cen A jSitarek fc Bednarek! |2010| ). Therefore, the cascade 
processes can be important in these sources provided that 
injection of primary particles occurs relatively close to the 
central engine as suggested by the recent observations of 
TeV 7-ray emission from some objects (e.g. Cen A, M87). 

We have already investigated the 7-ray production in 
blazars in a scenario of IC e pair cascade initiated by pri- 
mary 7-rays produced e.g., in the Synchrotron self-Compton 
model. These 7-rays develop IC e pair cascade in the ther- 
mal radiation fr om the accretion disk in t he whole volume 
above the disk (I Sitarek fc Bednarek! l2010l ). We have calcu- 
lated the 7-ray spectra produced in such a model for arbi- 



trary location of the observer in respect to the direction of 
the jet. Furthermore the influence of a magnetic field above 
the accretion disk on the development of the IC e pair 
cascade has been also investigated. 

In the present paper we assume that the cascade pro- 
cess in the accretion disk is initiated by leptons (electrons 
and positrons) acc elerated in the jet as postul ated by the 
external IC model (|Dermer fc Schlickeiser|[l993l ). These lep- 
tons produce the first generation of 7-rays by scattering an 
anisotropic disk radiation. Therefore the angular distribu- 
tion and spectrum of those primary 7-rays is dependent 
on the disk/jet system properties. They significantly differ 
from th e spectra obtained in term s of the previous calcu- 
lations jSitarek fc Bednarek! I2OI0I ) in which primary 7-rays 
were injected isotropically in the blob frame. These first gen- 
eration 7-rays initiate anisotropic IC pair cascade in the 
radiation field of the accretion disk. In a framework of this 
model we investigate spectral, angular and also time prop- 
erties of the 7-ray emission. Therefore, our present calcu- 
lations follow significantly more complicated and complete 
model than offered by the previously considered scenario 
with the isotropic injection of primary 7-rays produced in 
the SSC model. 



2 A MODEL FOR GAMMA-RAY 

PRODUCTION IN IC e ± PAIR CASCADE 

We adopt a classical optically thick and geometrically thin 
accretion disk model around super massive black hole as a 
dominant source of radi ation in the central region of the 
active galactic nuclei (see lShakura fc Sunvaevl (|l973l )). The 
emission of the accretion disk is treated as a black body with 
a power law temperature dependence on the distance, r, 
from the black hole, T = T in (r /r in )" 3/4 , where T in is the 
temperature at the inner disk radius ri n . We realize that the 
radiation field around the accretion disk can be much more 
complicated than produced in such a simple disk model. For 
example, the inner region of the accretion disk can become 
significantly hotter and geometrically thick as expected in 
some modern models for the advection dominated disks (e.g. 
iNaravan" fc Yil (|1994| )). A lso the disk may not be g eometri- 
cally thin, e.g. slim disks (|Abramowicz et al1ll988l ). 

This limitation of our calculations to such a simple disk 
scenario is motivated by the complexity of the overall cas- 
cade model which we intend to study in this paper. Other 
versions of the model, with modified accretion disk radiation 
field and/or the additional components of the radiation field 
produced e.g. isotropically around the disk, can be studied 
in the future along the prescription described in this paper. 

We assume that this disk thermal radiation create the 
only target for relativistic primary electrons. The blob is as- 
sumed to be point-like, i.e. its dimensions are much smaller 
than the characteristic distance scale (e.g. the injection 
height). The case of an extended injection blob with a com- 
plicated geometry can be also studied by a simply inte- 
gration of our results for the point-like injection over the 
whole injection region. Unfortunately, such process will re- 
quire to introduce additional free parameters, which describe 
the structure of the injection blob. We leave such more com- 
plicated studies (but probably also more realistic) for the 
future paper. Primary electrons are confined in a compact 
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blob moving along the jet (perpendicular to the disk plane) 
with a constant velocity v. Directions of injected electrons 
are isotropic in the blob frame. Note however, that for v = c, 
due to the relativistic beaming, distribution of electrons in 
the disk frame is highly anisotropic, with more electrons 
moving in the direction of the jet. In the disk reference frame 
the energy of the electron depends on its instantaneous di- 
rection. 

Primary electrons with specific spectrum interact with 
the disk radiation in the inverse Compton (IC) scattering 
process producing high energy 7-rays. These 7-rays prop- 
agate, in principle, at an arbitrary angle to the disk axis 
(following the instantaneous direction of the parent electron 
at the moment of the IC scattering). If the acceleration of 
electrons occurs relatively close to the accretion disk, then 
the optical depths for produced primary 7-rays can have 
substantial values at specific directions. 

Recently, we studied the optical depths for 7-rays prop- 
agating at an arbitrary directions in respect to the accretion 
disk for a few well known sources s uch a s 3C 273, 3C 279, 
Cen A jSitarek fc Bednarekl I2OO8L l20ld V It became clear 
that for the parameters characteristic for these objects a 
part of produced 7-rays can be absorbed in the disk radia- 
tion creating e* pairs. Thus, the first generation of pairs 
is created outside the jet in the whole volume above the 
accretion disk. These secondary e pairs propagate in the 
radiation field and the magnetic field above the accretion 
disk producing synchrotron radiation and the next genera- 
tion of 7-ray photons. As a result, IC e pair cascade develop 
in the volume above the accretion disk in which the addi- 
tional role is taken by synchrotron process of secondary e 
pairs. The details of such a model, in which the cascade is 
initiated by primary 7-rays produced in the jet (e.g. in the 
Synchrotro n self-Compton scenario) h as been recently con- 
sidered by ISitarek fc Bednarekl |201Cf ). It was shown that 
synchrotron energy losses of secondary cascade e ± pairs can 
be important only close to the disk surface where the dipole 
magnetic field is relatively strong. Farther out of the disk, 
the magnetic energy density drops faster than the disk radi- 
ation energy density. Therefore, we neglect the influence of 
the magnetic field on the IC e pair cascade process in the 
present studies. 

In this paper we consider another version of that sce- 
nario. The cascade process is initiated not by primary 7-rays 
but by primary electrons accelerated in the blob moving rel- 
ativistically along the jet. We assume that electrons are in- 
jected with a constant rate into the blob which moves along 
the jet. Electrons are trapped within the blob without any 
substantial escape. They cool in the accretion disk radia- 
tion field producing 7-rays which escape the blob into the 
whole volume above the disk. Other scenarios could be in 
principle also considered, e.g. continuous injection of elec- 
trons only at the base of the jet for certain period of time, 
or their injection into the extended blob. Unfortunately, the 
possible variety of the basic assumptions on the injection 
model of electrons can not be limited at present stage of our 
knowledge by the theory of the acceleration process of par- 
ticles in the AGNs. Therefore, we limit our calculations in 
the present paper to the simple possible scenario, i.e. con- 
stant injection rate of electrons along the jet. We develop 
a numerical code which follow the IC pair cascade in 
the 3D volume above the accretion disk. The method of the 



numerical simulation of the cascade scenario is described in 
AppendixfA] This code is also able to follow the time struc- 
ture of the high energy radiation which escapes towards the 
observer. All cascade particles (7-rays, e pairs) are tracked 
in this code up to the point when they escape far away from 
the disk (where the radiation field is already very weak) or 
their energies fall below assumed threshold equal to 10 GeV. 

As noted above, electrons are isotropic in the blob 
frame. They are injected with a power law spectrum which 
spectral index is fixed to a constant value independent on the 
distance from the base of the jet. As an example, we apply 
the differential spectral index equal to —2 for the spectrum 
extending from the minimum electron energy equal to 10 
GeV, up to a given maximum energy. These electrons cool 
on the IC process which efficiency depends on the distance 
from the accretion disk. So then, their equilibrium spectrum 
in the blob also depends on the distance from the base of the 
jet. We assume that electrons are trapped within the blob 
during its propagation through the accretion disk radiation 
field. We neglect the possible escape of electrons from the 
blob at this stage of calculations. It would require a detailed 
knowledge of the blob parameters (e.g. its dimensions, the 
magnetic field strength within the blob) to define the precise 
mechanism of escape of the electrons. In the future more re- 
alistic model, this effect should be taken into account. We 
expect that the energy dependent escape should steepen the 
equilibrium spectra of electrons within the blob which will 
reflect in obtained 7-ray spectra. On the other hand, escap- 
ing electrons will propagate in a complicated way around 
the accretion disk (depending on the magnetic field struc- 
ture) producing additional 7-ray photons. We are planning 
to study this effect in a future paper. 

The equilibrium spectrum of electrons for different loca- 
tions within the jet is determined numerically. Two models 
for the maximum energy of injected electrons are considered. 
In the first model, the maximum energy of electrons in the 
blob is fixed on 10 TeV. We apply such a high value in order 
to have possibility of investigation of cascade process also 
in the Klein-Nishina regime. In the second, more realistic 
model, the maximum energy is obtained by balancing the 
energy gains of electrons from the acceleration mechanism 
and their energy losses on the synchrotron process occurring 
in the magnetic field of the jet. We note that the magnetic 
field inside the jet is expected to be significantly stronger 
than in the remaining volume above the accretion disk due 
to its strong confinement by the jet. The energy gains of 
electrons from the acceleration mechanism in the blob can 
be parametrized by: 

ill) = W-Rl ~ 10 13 £B eVs" 1 , (1) 

V / acc 

where £ is the acceleration coefficient, 7?l is the Larmor ra- 
dius of electrons E is the energy of electron, and B is the 
magnetic field strength (in units of Gauss) at the accelera- 
tion region. The synchrotron energy loss rate is 

(ir) = | cct tC/b7c ~ 8.6 x 10- 4 B 2 7c 2 eVs- 1 , (2) 

V / syn 

where (7b is the energy density of magnetic field, <jt is the 
Thomson cross section, c is the velocity of light, and y e is the 
Lorentz factor of electrons. By comparing Eqs. [T] and [2l we 
get the limit on the maximum energy of accelerated electrons 
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due to the synchrotron losses: 



3.1 Blazars at large inclination angles (Cen A) 



E™* « 5(C- 2 /B) 1/2 TeV, 



(3) 



where £ = 0.01£-2. In a similar way we can obtain the limit 
on the maximum electron energy by comparing their energy 
gains with energy losses on IC process. The Klein-Nishina 
cross section drop at higher energies (KN effect). Therefore, 
unless the energy density of disk radiation do not dominate 
completely over the energy density of the magnetic field in 
the jet, this limitation is less restrictive. 

We estimate the magnetic field in the jet by assuming 
some level of the equipartition between the magnetic field 
and radiation at the disk inner radius, i.e Ub = Tj^iadj 
where 77 is a parameter which describes the level of equipar- 
tition. Based on this equality, we simply relate the mag- 
netic field at ri n to the inner disk temperature, Ti n by, 
Bin ~ 40^T 4 2 G, where T in = 10 4 T 4 K. If the blob is filling 
the cross section of a conical jet, we estimate the value of 
the magnetic field along the jet axis, 



B(h) = B in r in /(r in + h sin a) 



(4) 



The opening angle of the jet a is assumed of the order of 
~ 5°. Those magnetic fields are strong enough to confine 
primary electrons inside the blob. For such magnetic field 
model in the jet, we calculate the maximum energies of elec- 
trons at different locations along the jet (see Eq. [3|. In the 
cascade calculations we assume £ = 0.01. The first model 
for electron injection (model I) corresponds to the situation 
of a weak magnetic field in the jet (in respect to radiation 
field, i.e. rj « 1). Another possibility is the acceleration 
of particles in the rectilinear reconnection regions along the 
magnetic field lines without significant synchrotron energy 
losses. The second model (model II) corresponds to strong 
magnetic field in the jet. For the calculations performed in 
the framework of the second model we assume rj = 1. 



3 GAMMA-RAY SPECTRA EXPECTED AT 
DIFFERENT INCLINATION ANGLES 

The details of calculations of 7-ray spectra produced in 
such an anisotropic IC e pair cascade are described in 
Appendix [A] The injection spectrum of primary electrons, 
its cut-off energy at specific distances from the base of the 
jet, and other assumptions of the considered model are de- 
scribed in Sect. 2. The calculations are shown for the case 
of electron injection spectra normalized to 1 erg in the blob 
reference frame. The spectra of cascade 7-rays escaping to 
the observer will certainly depend on the observation an- 
gle (measured in respect to the jet axis) and the velocity 
of the blob. Therefore, we consider separately the cases of 
active galaxies with slower blobs but with the jets inclined 
at relatively large angles (e.g. such as Cen A, M87) and with 
faster blobs but inclined at relatively small angles to the line 
of sight. We note that the 7-ray cascading processes should 
occur efficiently, since the 7-spheres for these sources can 
extend to large distance s from the center of the ac cretion 
disk (see Figs. 4 and 5 inlSitarek fc Bednarek! (|2008h for 3C 
279 and Fig. 1 in lSitarek fc Bednarek! (|201Ch for Cen A). 



We apply the parameters expected for the misaligned blazar 
Cen A, as an example of an AGN seen at large inclina- 
tion angle. Cen A was observed in GeV 7-rays by the 
EGRET instrument on board of t he Compton GRP. T he 
measured spectral index of —2.58 (jHartman et alj 1 1999ft is 
rather soft. The emissio n was weak, with no significant flares 
ijSreekumar et ai]| 1999ft . Observations with the LAT instru- 
ment on board of the Fermi sa tellite show simila r behavior 
with a spectral index of —2.71 (|Abdo et al.ll2010l ). Cen A is 
the second radio galaxy (after M 87) observed in the VHE 
gamma-rays. With the integral flux of 0.8% of Crab Neb- 
ula flux, it is o ne of the weakest VHE gamma-ray sources 
observed so far (|Aharonian et al]|2009i ). The spectral index 
of VHE gamma-rays (—2.7) is consistent with the values 
obtained at the GeV energies. 

It is supposed that the central engine of Cen A har- 
bors a super massive black hole w ith the mass estimate d 
on McenA = (5.5 ± 3.0) x 10 7 M Q (ICappellari et al J 120091 ). 
A clear jet is observed in this radio galaxy. It propa- 
gates at the inclination a ngle estimated on aj e t ~ 15° — 
80° (|Horiuchi et all 120061 ). We apply the inner disk tem- 
perature equal to Ti n = 3 x 10 4 K, which is of the or- 
der of that one observed directly in other AGNs (e.g. 3C 
273). The total disk luminosity with these parameters is es- 
timated on Ld = 4'7rcrsB»'f n T i 4 1 1=3 3.6 x 10 41 erg s _1 , where 
r in = 4.5 x lO 5 A/ C cnA/M cm = 2.5 x 10 13 cm. Therefore, 
the disk luminosity is comparable to the measured nuclear 
X-ray luminosity p robably originating within the inner jet 
|Evans et al.112004 ). Note that the disk is observed at a rel- 
atively large angle, which effectively reduces the observed 
disk luminosity. This accretion disk emission can be addi- 
tionally obscured by the circumnuclear disk. Therefore, it 
may not be directly o bserved. The Cen A jet velocity is es- 
timat ed on P ~ 0.5c (|Tingav et all 1 19981 ; lHardcastle et al.l 
l2003h . 

In the case of blazars inclined at large angles (a > 10°), 
we consider blobs moving with mild velocities along the jet 
(i.e the Lorentz factor of the blob, 7b, is not much larger than 
1). As noted above, electrons in the blob frame are isotropic 
having a power law spectrum. As an example, we consider 
the differential spectral index equal to —2 above the min- 
imum electron energy equal to 10 GeV and the maximum 
energy described by the model I or model II (see Sect. 2). 
Injected electrons cool down on the IC process in the accre- 
tion disk radiation field when moving with the blob along 
the jet axis. The first generation of 7-rays initiate IC e pair 
cascades in the whole volume above the disk. 

We consider specific cases of the blobs propagating 
along the jet with the velocity v, starting at different height, 
Ho, above the accretion disk. At first, we consider a case 
of electrons originally injected in the blob at fixed distance 
Ho = 3, 10, 100ri n from its base. The blob is static or propa- 
gates along the jet with the velocity v = 0.5, 0.9c. Electrons 
cool on the IC process when propagating inside the blob. 
In Fig. Q] we show the 7-ray spectra escaping at different 
range of observation angles from such IC e pair cascade. 
Let us note the main features of these cascade 7-ray spectra. 
When electrons are injected close to the disk (see the case 
Ho = 3ri n ), the cut-offs in the 7-ray spectra appear at en- 
ergy for which the optical depth for 7-rays is comparable to 
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Figure 1. Gamma-ray spectra from the IC pair cascade ini- 
tiated by electrons injected with differential spectral index —2 
between 10 GeV and 10 TeV (in the blob frame). Electrons are 
injected isotropically in the blob starting at a distance from its 
base equal to: 3, 10, 100 ri„ (figures from left to right) with a 
velocity 0, 0.5c, and 0.9c. They cool gradually on the IC process. 
The 7-ray spectra are shown at the range of the observation an- 
gles § < 40° (solid), 40° < ^ < 60° (dashed), 60° < £ < 75° 
(dotted), and 75° < £ < 90° (dot-dashed). The calculations have 
been performed for the disk parameters as expected in the case 
of Cen A (see text). 



unity. This cut-off clearly depends on the observation angle 
£. It shifts to lower energies with increasing angle. This is 
due to a more efficient absorption effects, since 7-rays can 
interact at larger angles with more energetic soft photons 
coming from the bright, central part of the accretion disk. 
In the case of electrons injected farther from the black hole 
(e.g. Ho > 100), the 7-ray spectra extend clearly through 
the TeV energy range. The 7-ray spectra produced at larger 
inclination angles and the intermediate distances from the 
black hole show clear break. It corresponds to the transition 
between the Klein-Nishina and the Thomson regimes of IC 
scattering. In contrast, there is no break present for 7-ray 
spectra in the case of a small £ and large Ho- In this case 
electrons are not able to interact efficiently with the higher 
energy photons coming from the inner part of the accretion 
disk. The final cut-offs in the 7-ray spectra are related to the 
maximum energies of primary electrons in the disk frame. 

For blobs moving with small velocities, the GeV 7-ray 
flux emitted at larger angles can even dominate over the 
7-ray flux emitted at lower angles. This is caused by more 
efficient cascading process for photons which propagate at 
large inclination angles to the disk axis. For faster jets, the 
effects of relativistic beaming start to become clearly visi- 
ble. In this case, directions of electrons in the disk reference 
frame are more preferably pointed along the jet. Therefore, 
also the 7-rays are produced more preferably along the jet 
direction. Due to this beaming effect, the 7-ray flux observed 
at higher inclination angles is much lower than that one ob- 
served along the jet. Moreover, energies of 7-ray photons 





-v = 0.5c, H = 


1 - 100r in - 


2 


- v = 


0.9c, H = 1 - 100r in - 


. *• , 


s 

• M 

- Ni ~ 


1 


-1 




•* \ v \ - 

* **. N 1- 

s '•. > 1: 
• "'*• 1 l : 
^ *'•• 1 Y 




1 ': 1 ; 


-2 








»v\ 


-3 




• l f 
\\l ] 

., 1 ,.,,!,... T ,,,, 1 .'. ,,. 1 


1.5 2 2.5 3 


3.5 4 4.5 




1.5 


2 2.5 3 3.5 4 4.5 


- 7 v = 0.5c, H = 


1 - 300r in -_ 


2 


r v = 


0.9c, H = 1 - 300r in -_ 


is - s ■ m • ^, >r ^ 




1 




""*••» "V \ - 


- s, " - 













•""••.\\ 1 
s '■ 1 

i \ 


-1 




* - 4 1 
x i » \- 






-2 








t 1 : 


-3 




H » \ 



1 1.5 2 2.5 3 3.5 4 4.5 1 1.5 2 2.5 3 3.5 4 4.5 
log (E GeV) log (E/GeV) 

Figure 2. As in Fig. [T] but for the extended injection region of 
relativistic electrons along the jet. It is assumed that the injection 
process occurs between l-30ri n (top figures), l-100ri n (middle) 
and l-300ri„ (bottom). The blob propagates with the velocity: 
0.5c (left figures) and 0.9c (right figures). The 7-ray spectra are 
shown within the range of observation angles: £ < 40° (solid), 
40° < £ < 60° (dashed), 60° < £ < 75° (dotted), 75° < £ < 90° 
(dot-dashed). 



moving along the jet (observed in the disk reference frame) 
are relativistically boosted. As a result, the cut-offs of cas- 
cade 7-ray spectra appear at higher energies. 

In Fig. [2] we show the 7-ray spectra expected at differ- 
ent observation angles assuming that primary electrons are 
injected into the blob when it moves through a range of dis- 
tances above the accretion disk between Hi and H2- Their 
cooling process is followed up to distances at which they ob- 
tain energies below 10 GeV. The 7-ray spectra are shown for 
v = 0.5, and 0.9c, H ± = lr in and H 2 = 30, 100, 300r in . As 
expected from the analysis of the spectra calculated for fixed 
injection place (see Fig. [l|, the shape of the 7-ray spectra 
strongly depends on the range of injection distances of elec- 
trons. The 7-ray spectra extend to higher energies for the in- 
jection of primary electrons farther from the accretion disk. 
They are also flatter. Still, the GeV 7-ray fluxes produced 
at large inclination angles dominate for slow jets. However, 
the 7-ray spectra at the TeV energies become steeper. They 
also have a lower energy cut-off for large inclination angles. 
Slower jets produce lower 7-ray fluxes with steeper spectra 
at TeV energies due to a weaker relativistic beaming and a 
more significant KN effect. We conclude that the location of 
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Figure 3. Comparison of cascade gamma-ray spectra (thin 
curves) with the observations of Cen A. The IC e^ pair cas- 
cade is initiated by electrons injected with a power law with an 
index —2.5 (left figures) and —3 (right) inside the blob propagat- 
ing with a velocity: 0.5c. The injection process occurs between 
l-100ri„ (top figures), and l-300ri n (bottom). The 7-ray spec- 
tra are shown within the range of observation angles: £ < 40° 
(thin solid), 40° < $ < 60° (thin dashed), 60° < § < 75° (thin 
dotted), 75° < § < 90° (thin dot-dashed). Thick curves (with 
error bow-ties) show th e observations of Cen A: EGRET (d otted, 
ISreekumar et alJll99Sh. Fermi (dashed . lAbdo et alTl2010h . and 
H.E.S.S. (solid. lAharonian et alj|2009h The cascade spectra are 
normalized to the total energy output in the 0.5 - 5TeV range as 
measured by the H.E.S.S. telescopes. 
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Figure 4. As in Fig. [2] but for the model II of the maximum 
energy of electrons injected into the blob. These maximum ener- 
gies depend on the distance along the jet according to Eq. \3\ The 
magnetic field inside the jet is defined by r\ = 1. 



the break in the 7-ray spectrum and its spectral index above 
~ 100 GeV is clearly related to the inclination angle of the 
jet. It may serve as a diagnostic tool for the inclination of 
the accretion disk with respect to the observer. On the other 
hand, the information on the inclination angle can allow to 
estimate a possible break in the expected 7-ray spectrum 
from a specific source. 

In Fig. Owe compare the calculated cascade 7-ray spec- 
tra with the observations of Cen A. The velocity of the blob 
is fixed on v = 0.5c (as in lHardcastle et all l|2003h ). Electron 
are injected with the spectral index —2.5, and —3. Reason- 
able agreement with the spectral shapes observed by the 
EGRET and the H.E.S.S. is obtained for low observation 
angles (< 40°) and the spectral index of injected electrons 
equal to —3. The 7-ray spectra produced in our model for the 
larger inclination angles of the accretion disk-jet system are 
too steep to be able to explain these measurements. Note 
however, that the results reported by the EGRET, Fermi 
and H.E.S.S. are not simultaneous. So then, this conclusion 
should not be considered as a final. 

The 7-ray spectra obtained for electrons with maximum 
energies limited by synchrotron energy losses (model II) are 
shown in Fig. [4] In order to consider more extreme example 
of model II, we apply the synchrotron energy losses calcu- 
lated for the value of r\ — 1. For this value of 77, the cut-off 
in the primary electron spectrum is at lower energy than in 



the case of model I. It results in slightly steeper 7-ray spec- 
tra with the cut-offs at lower energies. The effect becomes 
more pronounced in the case of a blob propagating closer to 
the black hole, due to stronger magnetic fields and therefore 
lower maximum energies of accelerated electrons. 



3.2 Blazars at small inclination angles (3C 279) 

As an example of a blazar observed at a small inclina- 
tion angle we take the parameters of the famous blazar 3C 
279. With a redshift of z = 0.538, it is a rather distant 
sourc e. It is one of the ea rliest discovered blazar in the 7- 
rays (|Hartman et al.lll992l ) . A strong 7-ray flare with a very 
flat spectrum (differential spectral index —1.89) has been de- 
tected from 3C 279 already in June 1991. The flux increased 
during ~ 1 week period an d declined during the following 
2 days (|Kniffen et al.lll993l ). Even stronger 7-ray flare has 
been observed by the EGRET in 1996 with a sim ilar time 
structure and spectral index (|Wehrle et al.l 11998) ^ . 3C 279 
has been positively det ected by the EGRET telescope in 
all observation periods (|Hartman et al.l l2001al ). showing a 
variety of flux and spectral index stages. The shortest vari- 
abili ty time scales obse r ved in these flares wer e below ~ 1 
day (|Wehrle et al.lll998l ; iHaxtman et alj|2001bl ). 

More recent observations performed with the LAT tele- 
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Figure 5. Gamma-ray spectra from IC e pair cascade initiated 
by electrons injected isotropically in a blob with spectral index 
—2 between 10 GeV and 10 TeV (in the blob frame) for disk pa- 
rameters as expected for 3C 279. A blob propagates along the 
jet starting from a distance 3, 10, and 100r; n from its base (pan- 
els from left to right), with a velocity 0.9c (top panels), 0.99c 
(middle), and 0.998c (bottom). The spectra are shown for the 
range of inclination angles £ < 10° (solid), 10° < § < 20° 
(dashed), 20° < £ < 30° (dotted), 30° < £ < 40° (dot-dashed), 
40° < £ < 50° (dot-dot-dashed), and § > 50° (dot-dot-dot- 
dashed). 



scope (high energy gamma-ray instrument on board of the 
Fermi) show a different behavior. In December 2008 a longer 
flare, with a ~ 15 days period of constantly rising flux fol- 
lowed by a slower decay lasting for ~ 30 days, was reported 
dur ing period of a high activ ity (the LAT monitored source 
list. lCiprini fc Chatv I l|200Sl )). A similar, long flare has been 



obser ved by Fermi-LAT in July /August 2009 (jlafrate et al.l 
l2009h . 

3C 27 9 has been also rece ntly detected by the MAGIC 
telescope ((Albert et al.l l2008bl ). being at present the most 
distant sub- TeV 7-ray source. Due to a severe absorption 
in the extragalactic background light (EBL) radiation, the 
observed spectrum in sub- TeV range is very soft, with a 
spectral index of —4.11 ± 0.68. 

For this source we apply the disk luminosity Ldisk = 
2 x 10 45 ergs _1 a nd the temperatu re Tin = 2 x 10 4 K at the 
disk inner radius (IPian et alj|l999h . The inner radius of the 
accretion disk in 3C 279 is estimated on ri n = 4.2 x 10 15 cm. 



3.2.1 Spectra at the source 

As in the case of Cen A, let's at first consider a simple sce- 
nario of the injection source of electrons at fixed distance 
from the base of the jet. We assume the injection spectrum 
of electrons has a differential spectral index —2 between 
10 GeV and 10 TeV. In Fig. [5] we present the angle de- 
pendent 7-ray spectra produced by electrons injected at a 
distance from the accretion disk Ho = 3, 10, 100 ri n . The 
calculations have been performed for various values of the 
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Figure 6. As in Fig. [5] but for the extended injection region 
of relativistic electrons along the jet. The injection of electrons 
occurs between l-30ri„ (top figures), l-100r; n (middle), and 1- 
300ri„ (bottom). The blob propagates with a velocity 0.99c (7 = 
7.1, left figures) and 0.998c (7 = 15.8, right figures). 



blob velocities v = 0.9c (corresponding to the Lorentz fac- 
tor 7 = 2.3), 0.99c (7 = 7.1), and 0.998c (7 = 15.8). These 
Lorentz factors are more characteristic for distant blazars 
observed at small inclination angles. Those sources can be 
detected due to the strong relativistic beaming, which signif- 
icantly enhances the observed flux. The shape of the cascade 
7-ray spectra remains rather stable, with a cut-off energy 
determined by the injection height and the observation an- 
gle. We note, that the change of the spectral index due to 
the KN effect is not so clearly observed as in the case of 
Cen A. It is only visible in the 7-ray spectrum produced by 
electrons injected at the distance 100ri n from the disk. This 
much smaller KN effect is due to the lower energies of the 
disk photons seen by electrons in the reference frame of the 
blob moving with significantly larger Lorentz factors than 
in the case of Cen A. 

We also calculate the 7-ray spectra in the case of a con- 
tinuous injection of electrons into the blob moving through 
a range of distances within the jet. The 7-ray spectra es- 
caping in this case at different range of angles are shown 
in Fig. [6] The 7-ray spectra extend to lower energies for 
larger inclination angles due to stronger absorption effects 
and also kinematic effects caused by lower Doppler factors 
(D = [7b ( 1 — /9b cos £)] ~ 1 ) . For fast blobs, the 7-ray emission 
is strongly limited to the small inclination angles (£ < 10°). 
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The cut-offs in the 7-ray spectra also clearly depend on 
the location of the injection region of electrons within the 
jet. However, the spectral index in the lower energy part of 
the 7-ray spectra (the GeV energy range) does not depend 
strongly on the observation angle. 

3. 2. 2 Spectra after propagation in the EBL 

Since 3C 279 is a distant source, the observed 7-ray spectra 
differ from the spectra at the source presented in Fig. [6] Due 
to the cosmological effects, energies of 7-ray photons are 
reduced. Moreover, the part of the 7-ray spectrum at a few 
hundred GeV is strongly absorbed in the EBL. In order to 
obtain spectra and light curves as seen by a distant observer 
one has to correct for these propagation effects on the way 
from the source to the observer. The corrections depend on 
the redshift 2 of the source. The observed energy E b s of the 
7-ray photon is reduced according to E \^ s — E cm /(1 + 2), 
where _E em is the energy of the emitted 7-ray in the disk 
frame. The observed 7-ray spectra and the light curves are 
attenuated by a survival probability e~ TEBL , where tebl is 
the energy dependent optical depth in the EBL radiation. 
The cosmological effects also influence the arrival time of 7- 
ray flare according to At b s = At cm (l + z), where A em is the 
time interval measured in the disk frame and At b s is the 
time interval measured in the observer frame. We took all 
those effects into account when obtaining the 7-ray spectra 
in Fig. □ 

The 7-ray optical depths in the EBL radiation field , 
tebl, are calculated according to iFranceschini et all i|2008l ) 
model. Due to a large redshift of 3C 279, the 7-ray spectra 
above ~ 200 GeV are strongly absorbed resulting in a strong 
steepening of the spectrum. Therefore, the shape of the 7- 
ray spectra above ~100 GeV are quite independent on the 
blob Lorentz factor and the extend of the injection region 
within the jet. This creates problems for extracting physical 
parameters of the 7-ray emission region from the comparison 
of the observed and calculated spectra. 

The IC e pair cascade 7-ray spectra, after absorption 
effects in the EBL, are compared with the observations of 
the two high states of the source seen by the EGRET and 
the MAGIC, and with mean flux level spectrum seen by the 
Fermi-L AT (see Fig. [8]). The spectrum calculated for the 
beginning of the flare (t = — 23 days) and for the range 
of injection distances Hq — l-300ri„ can fit the EGRET 
and the MAGIC measurements provided that electron are 
injected with a power law spectrum and spectral index in the 
range from —2.5 to —3. Note that the model predicts how 
the level of emission should evaluate in time. After ~ (60-90) 
days, it should drop to the mean emission level determined 
by the Fermi-LAT observations. Thus a single large flare 
every several months can explain the mean emission level. 



4 TIME STRUCTURE OF THE GAMMA-RAY 
EMISSION 

As shown in the previous section the disk radiation field seen 
by electrons in the blob moving along the jet changes signif- 
icantly. As a result, 7-rays are produced at different places 
above the disk and at different moments. Therefore, the 7- 
ray spectra arriving to the distant observer show significant 




Figure 7. Gamma-ray spectra from 3C 279 for the parameters as 
in Fig. [6] but after the prop agation in extragalactic b ackground 
light described by the model [Sanccschini et al. (200j|). 




Figure 8. The IC pair cascade 7-ray spectra at specific time 
intervals calculated for 3C 279 after the propagation in the inter- 
galactic radiation, compared with various observations. Primary 
electrons are injected with the power law spectrum and spectral 
index —2.5 (left figure) or —3 (right figure) between 10 GeV and 
10 TcV in the blob frame of reference. The injection occurs in the 
range of distances from the base of the jet: Hq = 1— 300ri n . 7-ray 
spectra are shown for the time intervals: t = 0—23 days (solid), 
23-45 days (dashed), 45-68 days (dotted), 68-90 days (dotted). 
Blob propagates with a velocity v = 0.998c (7 = 15.8). Thin lines 
with error bow-ties show observational re sults in a high stat e 
of 3C 279 collected by the MAGIC (solid. lAlbert et al.ll2008'bT) . 
the EGRET (black points, lHartman et al ]|2001at) , and the mean 
spectrum from 11 month observations by the Fermi-LAT (dashed, 
lAbdo et al .|2010T) . The solid curve is normalized to the 7-ray spec- 
trum observed by the MAGIC telescope. 
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dependence on the observation time. Here we calculate the 
7-ray light curves expected from such IC e pair cascade 
for specific parameters of the model and as a function of the 
location of the observer in respect to the disk axis. 7-rays 
are produced in a cascade initiated by electrons which are 
injected with a power law spectrum continuously along the 
jet and cool in the IC process (see for details Section [3]). 
As before, we consider separately the cases of slowly moving 
jets, but observed at a large inclination angle, and fast mov- 
ing jets, observed nearly along the jet direction. For each 
case we calculate the expected 7-ray light curves and spec- 
tra arriving to the observer at different moments for specific 
propagation distances of the blob with fixed velocity. The 7- 
ray light curves calculated for different parameters are nor- 
malized to I erg of injected energy of electrons (in the blob 
frame) . 

4.1 Blazars at large inclination angles 

In order to investigate the time dependence of 7-ray emission 
let's first consider a simple scenario in which electrons are in- 
jected at a fixed distance from the accretion disk. Electrons 
with a differential power law spectrum are instantaneously 
injected inside a compact blob at a given height Ha above 
the accretion disk. The blob is either at rest, or it is mov- 
ing along the jet with a velocity v. The electrons are cooled 
down in the accretion disk radiation field producing 7-rays, 
which initiate the IC e pair cascade in the whole volume 
above the accretion disk. These 7-rays arrive to the observer 
at different times from two reasons. At first, primary 7-rays 
are produced at different distances from the base of the jet. 
Afterwards, secondary cascade 7-rays are produced at dif- 
ferent locations above the accretion disk passing different 
distances on their way to the observer. 

In Fig. [9] we show the 7-ray light curves obtained in 
such a scenario for the blob at rest or moving with the ve- 
locity v = 0.9c, and for two injection heights along the jet 
Ho — 10 and I00n n . In most of the cases, the 7-ray light 
curves show a strong initial peak with an exponential de- 
cay. The characteristic decay time strongly depends on the 
injection height. If the injection place is located far away 
from the disk, electrons cool down in a relatively weak ra- 
diation field. Therefore, cooling process of electrons takes 
longer producing broader peak (see right figures in Fig. [9]) . 

It is interesting to evaluate the effect of the blob velocity 
on the 7-ray light curves. For an injection place of electrons 
close to the disk, the emission peak is sharper in the case of a 
faster blob and low inclination angles, i.e. the flare is charac- 
terized by shorter time scales. This effect is connected with 
the relativistic beaming of the primary 7-rays along the jet 
direction. During typical cooling time scale of electrons, the 
blob moves in the direction of the observer, effectively re- 
ducing the time scale in the reference frame of the observer. 
In the case of the electron injection place located farther 
away from the disk, the 7-ray light curves show a long tail 
for a fast blob. This is caused by the comparable time scales 
for the electron cooling in the disk radiation with the travel 
time scale of the blob along the jet. Therefore, the 7-rays are 
produced on much longer distances and the corresponding 
7-ray emission time scale in the observer's reference frame 
increases. It results in a long tail with a slower than expo- 
nential decay. 
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Figure 9. The 7-ray light curves from the IC pair cascade 
initiated by electrons, injected isotropically in a blob, with spec- 
tral index —2 between 10 GeV and 10 TeV (in the blob frame) for 
Cen A at different energy ranges E > 10 (solid), > 30 (dashed), 
> 100 (dotted), and > 300 GeV (dot-dashed). Injection of elec- 
trons occurs at H = 10 r; n (left figures), and 100 r- ln (right fig- 
ures). Stationary blob (v = 0) (top figures) and blob propagating 
with v = 0.9c (bottom figures). Two observation angles are con- 
sidered 32° (thick lines) and 57° (thin lines). 



In the next step, we consider the injection of the elec- 
trons in the blob at a fixed range of distances from the 
disk. We are considering both the cases of a slow (velocity 
v — 0.5c) and a fast (v = 0.9c) moving blob. It is assumed 
that the injection process of electrons occurs in the range of 
distances between lr m and 30, 100, and 300 n n . 

The results of calculations are shown in Fig. 1101 In gen- 
eral the 7-ray light curves show at first a plateau with a 
strong exponential decay. Depending on the injection pa- 
rameters, it can be followed by a long tail of slowly decaying 
emission. The flat part of the light curve is produced by the 
fast cooling electrons close to their injection place. Depend- 
ing on the range of injection heights, the emission can last 
at a fairly constant level for from a few up to a few tens of 
hours. The tails, as in the case of a fixed injection, are caused 
by electrons escaping into region of a constantly decreasing 
radiation field. The duration of the flat part of the light 
curve is determined mostly by the range of the distances on 
which the injection of electrons takes place. Moreover, for a 
fast moving blob observed at low inclination angles, it is sig- 
nificantly reduced due to the relativistic effects. If a source 
is observed at large inclination angles, then the relativistic 
effects are less pronounced, resulting in a weaker but more 
extended plateau in the 7-ray light curve (see thin curves in 
right figures in Fig. [T0)l . 

The properties of the light curve depend on the en- 
ergies of 7-ray photons. The beginning of the flare, corre- 
sponds to the time in which the blob is close to the disk. In 
such a strong radiation field, the very high energy (VHE) 7- 
rays are strongly attenuated. Therefore, the VHE emission 
(> 300 GeV) is weaker at the beginning of the flare (see 
dot-dashed curves in top Fig. 1 10 p . The absorbed 7-rays are 
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Figure 10. The 7-ray light curves produced in IC pair cas- 
cade initiated by electrons injected isotropically in a blob. Elec- 
trons are injected with a spectral index —2 between 10 GeV and 
10 TeV (in the blob frame). The parameters of Cen A are ap- 
plied. The injection process of electrons occurs between 1-30 rj n 
(top figures), 1-100 rj n (middle figures) and 1-300 ri„ (bottom 
figures). The blob propagates with the velocity: 0.5c (left figures) 
and 0.9c (right figures). The light curves are shown for energy 
thresholds and angles as in Fig. [9] 



reprocessed in the IC e pair cascade into a lower energy 
7-rays. Those secondary 7-rays create an additional peak in 
the light curve above 10 GeV at the beginning of the flare 
(see top left Fig. I10[l . This leads to an interesting conse- 
quence. VHE 7-ray emission (in the Cherenkov telescopes 
energy range) should appear after the GeV 7-ray emission 
(detected by 7-ray satellites). Therefore, the model predicts 
a soft to hard 7-ray evolution in the case of 7-ray flares 
observed from blazars. This effect is even more pronounced 
for large observation angles (see dot-dashed thin lines in 
Fig. UOp . In this case, due to a more favorable angle be- 
tween the direction of the 7-ray photon and the hot inner 
region on the accretion disk, the absorption of high energy 7- 
rays is even stronger. Thus, primary TeV and sub- TeV 7-ray 
photons are reprocessed more efficiently into the multi-GeV 
range. 

In Fig. [TT]we show the 7-ray spectra for different time 
intervals in the case of the model I for the maximum en- 
ergy of electrons. In the case of electron injection close to 
the disk, the multi-GeV 7-ray flux is slowly dropping with 
time in the plateau phase. However, the spectrum extends 
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Figure 11. The 7-ray spectra at a specific time intervals t from 
the IC pair cascade initiated by electrons, injected isotropi- 
cally in a blob, with a differential spectrum with a spectral index 
—2 between 10 GeV and 10 TeV (in the blob frame) for the pa- 
rameters of Cen A. Blob propagates with the velocity 0.5c (left fig- 
ures) and 0.9c (right figures). Observation angle 32° (thick lines) 
and 57° (thin lines). The 7-ray spectra are accumulated in the 
time intervals: t = 0-2h (solid curves), 2-4h (dashed), and 4-6h 
(dotted) for the injection heights in the range Ho = l-30ri n (up- 
per panel); t = 0-10h (solid), 10-20h (dashed), 20-30h (dotted) 
for H = l-100r;„ (middle panel); and t = 0-50h (solid), 50-100h 
(dashed), 100-150h (dotted) for H = l-300r iri (bottom panel). 



to higher energies with the development of the flare reaching 
TeV energies after a few hours (see the upper left Fig. 
On the other hand, for the injection region of electrons ex- 
tending far away from the disk, the 7-ray spectrum keeps 
a constant shape in time. The level of emission decreases 
significantly while entering the tail phase in the time struc- 
ture of the emission (see dotted curves in bottom panels in 
Fig. Hip . This is a result of a slower cooling of electrons in 
the constantly weakening radiation field without strong ab- 
sorption of 7-rays produced in the IC process. We note that 
the higher value of blob velocity result in larger variations 
of the spectrum. 

In Fig. [T2] and Fig. 1131 we present the 7-ray light curves 
and the time dependent gamma-ray spectra for the maxi- 
mum energies of electrons obtained in terms of the model 
II. The basic properties of the 7-ray light curves and spec- 
tra are similar to those obtained in term of the model I. 
In general, electrons are accelerated to lower energies in the 
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Figure 12. As in Fig, 1101 but for the maximum energies of elec- 
trons obtained according to model II (see text). 
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Figure 13. As in Fig. 1111 but for model II of maximum energy 
of injected electron. 



model II. Therefore the spectra also extent to lower energies. 
The differences are most pronounced in the case of electron 
injection close to the disk. Then, the maximum energies of 
electron determined by the synchrotron energy losses are 
clearly below 10 TeV. This results in the 7-ray spectra lim- 
ited to lower energies than calculated in model I (see upper 
plots in Fig. El). 

4.2 Blazars at small inclination angles 

The 7-ray light curves from blazars observed at low inclina- 
tion angles are analyzed taking as an example the param- 
eters of 3C 279 (see Fig. I14p . In general, the light curves 
show at the beginning slow dependence on time (a plateau) 
followed by a strong exponential cut-off. In the case of faster 
blobs the cooling process of electrons occurs in a more ex- 
tended region along the jet. Then, the 7-ray light curves 
show longer and slowly decaying tails. If the injection re- 
gion is limited to the inner part of the jet, then the light 
curves show only a flat plateau extending for a few days. 

The duration of the 7-ray flare depends on the observa- 
tion angle. For small angles, the tail emission can even last 
a few times longer than the duration of the plateau (see e.g. 
the 7-ray light curves in the middle and bottom panel in 
Fig. I14[) . However, since the tail emission is dropping fast, 
depending on the total flare luminosity and the sensitivity 
of the instrument, only it's beginning may be observable. 



On the other hand, for larger inclination angles the 7-ray 
emission is much more confined in time. The 7-ray flares 
observed at lower energies lasts significantly longer than at 
higher energies. Therefore, we should expect shorter flares 
at TeV energies than at GeV energies. 

The 7-ray spectra at the source, expected at specific 
time intervals, are shown in Fig. 1151 In the case of elec- 
tron injection close to the accretion disk (injection range 
1-30 ri n ) the 7-ray spectra have similar shape during the 
plateau phase (see upper panel in Fig. I15|) . For electrons in- 
jected along larger range of distances from the base of the 
jet, the 7-ray spectra during the plateau and exponential 
decay phase differ significantly. During the plateau phase, 
the spectra has similar shape and slightly extend to higher 
energies with development of the flare. However, during the 
exponential decay phase, the 7-ray spectra become signif- 
icantly steeper with a cut-off at lower energies. We note 
that the 7-ray emission during the tail phase is produced by 
electrons that are in the pure cooling phase above the injec- 
tion region limited by the distance Hi- Therefore, spectrum 
of these electrons is clearly steeper than the assumed one, 
oc E~ 2 , in the acceleration region. 

The 7-ray light curves and spectra at the specific time 
intervals, after correction for the propagation effects in the 
intergalactic radiation field, are shown in Fig. [inland Fig. 1 171 
The time intervals in the Fig. [17] correspond to those shown 
in Fig. 1151 calculated for 3C 279 in its reference frame. The 
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Figure 14. The 7-ray light curves calculated for the electron 
spectrum and their injection distances as in Fig. 1101 but for the 
parameters of 3C 279. The blob is moving with the velocity v = 
0.99c (7 = 7.1, left figures), or 0.998c (7 = 15.8, right). The line 
styles correspond to these same energy ranges as in Fig. 1101 The 
values of the observation angles have been fixed on 12° (thick 
lines) and 22° (thin lines). 

clear differences between these spectra are caused by the 
cosmological propagation effects. Because of these effects, 
the observed time scales of the flares are also diluted by a 
factor of 1 + z with respect to the time measured at the 
source. 

The propagation in the EBL over the distance scale 
of 3C 279 results in an almost complete absorption of the 
7-ray spectrum above ~ 1 TeV and a very steep spectrum 
above ~ 100 GeV. The shapes of the 7-ray spectra at specific 
time intervals become more similar to each other showing 
weaker dependence on time with the development of the 
flare. Still, the spectral index at 100 GeV of the tail emission 
is significantly steeper than in the plateau emission phase. 

For the illustration, we confront the calculated 7-ray 
light curve for specific parameters with the recent observa- 
tions of 3C 279. The Fermi-LAT instrument revealed a large 
flare lasting for several weeks. The LAT 7-ray light curve 
above the energy threshold of 0.1 GeV r esulting from the 
LAT bright source monitoring (see also ICiprini fc Chatvl 
120081 ) has been used. This light curve has a similar shape 
as the energy light curve observed above 1 GeV, but due to 
larger statistics, it has a more precise sampling. It is com- 
pared with the 7-ray light curve at the energy threshold 



Figure 15. Time dependent spectra for 3C 279. Blob propagates 
with the velocity 0.99c (7 = 7.1, left figures) and 0.998c (7 = 15.8, 
right). The upper figures are for the electrons injected within 
the range of distances from the base of the jet Ho = l-30r;„ 
(spectra for times (in the source frame of reference) t = 0-1 day 
(solid curves), 1-2 days (dashed), 2—3 days (dotted), 3—4 days 
(dot-dashed)). The middle figures arc for Ho = l-100r; n (spectra 
for times t = 0—5 days (solid), 5—10 days (dashed), 10-15 days 
(dotted), 15—20 days (dot-dashed)). The bottom figures are for 
Hq = l-300r; n (spectra for times t = 0-15 days (solid), 15—30 
days (dashed), 30-45 days (dotted), 45—60 days (dot-dashed)). 
The observation angle is fixed on 12° (thick lines) and 22° (thin 
lines) . 



> 10 GeV calculated in the framework of our model (see 
Fig. ll8[) . The behavior of the dropping part of the 7-ray light 
curve from 3C 279 is consistent with our calculations for the 
observation angle equal to 9°. However, the rising part of the 
calculated light curve is very steep which is due to our as- 
sumption on the injection of electrons into a point-like blob 
at fixed moment of time. In the case of an extended blob or 
the injection mechanism operating at the base of the jet for 
a specific period of time with rising efficiency, the finite rise 
time of the 7-ray flare from 3C 279 should be naturally ex- 
plained in terms of our model. Therefore, we conclude that 
the several weeks 7-ray flares, observed in the GeV range 
from the blazars of the 3C 279 type, can be explained by 
our IC e pair cascade scenario initiated by electrons in the 
radiation field of the accretion disk. 
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Figure 16. The 7-ray light curves calculated for 3C 279 (for the 
parameters as shown in Fig. ll4| l but after the propagation through 
the intergalactic radiatio n field which is taken from the model by 
Franccschi ni et al] l|2008l ). The light curves are presented in the 
observer's frame on the Earth. 



5 DISCUSSION AND CONCLUSIONS 

We discuss the most complete version of the external IC 
model for the 7-ray production in blazars. It is assumed that 
electrons are injected in the inner part of the jet and comp- 
tonize thermal radiation coming from the accretion disk sur- 
face. As a result a first generation of 7-rays is produced. 
These 7-rays initiate IC e pair cascade in the anisotropic 
disk radiation. We calculate the 7-ray spectra escaping to 
the observer located at an arbitrary angle in respect to the 
direction of the jet. In the framework of this model we cal- 
culate also the 7-ray light curves expected from blazars 
in which blobs containing relativistic electrons propagate 
with arbitrary velocities. The angular dependent 7-ray light 
curves and spectra are calculated for blazars observed at 
large and small inclination angles. As an example we have 
applied the parameters of the two well known AGNs: Cen A 
and 3C 279, which have been recently observed in the TeV 
7-rays. The first source is located relatively nearby, and has 
a jet which propagates at a relatively large angle in respect 
to the observer with a mild velocity. In contrast, the second 
source is more distant. It's inclined at a small angle and the 
jet propagates with a large Doppler factor. 

In order to perform the cascade calculations with a lim- 
ited number of free parameters, we made some simplified 
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Figure 17. The 7-ray spectra at specific observed time intervals 
calculated for 3C 279 (for the parameters as shown in Fig. I15|l 
but after propagation in the intergalactic radiation. The observed 
times are 1 + z larger than the local times i.e.: top figures (the 
injection height range Ho = 1— 30ri n ) (spectra for times t = 0-1.5 
days (solid curves), 1.5-3 days (dashed), 3-4.5 days (dotted), 4.5- 
6 days (dot-dashed)). Middle figures H = 1-lOOr^ (t = 0-7.5 
days (solid), 7.5-15 days (dashed), 15—23 days (dotted), 23—30 
days (dot-dashed)). Bottom figures Ho = l-300ri„ (t = 0—23 
days (solid), 23-45 days (dashed), 45-68 days (dotted), 68-90 
days (dot-dashed)). 



assumptions. At first, we assume that the disk radiation 
is the only source of soft photons with which relativistic 
electrons can interact. The calculations with other radiation 
fields (e.g. infrared radiation from the molecular torus or the 
disk and jet radiation scattered by the matter around the 
accretion disk) are straightforward. But the interpretation 
of the final results will be more difficult due to the large 
number of free parameters describing such more realistic 
model. Therefore, the analysis of such an anisotropic cas- 
cade model with a more complicated soft radiation field is 
left for a future paper. We also assume the simplest possible 
(i.e. constant along the jet) dependence of the injection rate 
of relativistic electrons along the jet. This allows us to study 
the basic features of this complicated cascade model. These 
features may change significantly for the more complicated 
forms of the injection rates along the jet. In fact such more 
complicated dependence for the injection rate is suggested 
by the observed 7-ray light curves of flares observed from 
blazars at different energy ranges. The simultaneous future 
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Figure 18. The cascade 7-ray light curve (above 10 GeV) in 
the observer's reference frame (as in Fig. I16|l is confronted with 
the shape of the 7-ray flare observed from 3C 279 by the Fermi- 
LAT at energies above O.lGeV (the LAT monitored source list, 
ICiprini fc Chatv I fcOOSft l. Primary electrons are injected with the 
power law spectrum and spectral index —2.5 at a constant rate 
inside the blob moving along the jet with the velocity 0.998c 
through the range of distances from the base of the jet Ho = 1— 
300r; n . The observation angle is assumed to be 9° (solid curve) 
or 15° (dashed). 



observations of 7-ray light curves at different energy ranges, 
e.g. GeV and TeV, will put constraints on the detailed mod- 
els of 7-ray emission from specific blazars. 

We considered two different models for the maximum 
energies of injected electrons. In the model I, the maximum 
energies of electrons are fixed on 10 TeV independent on the 
distance from the base of the jet. In the second model, elec- 
tron maximum energies are determined by their synchrotron 
energy losses in the local magnetic field of the jet. The elec- 
trons have clearly lower energies in the second model. 

Based on our calculations, we predict a few character- 
istic features of the IC pair cascade model which can 
be tested by the future multiwavelength GeV and TeV 7- 
ray observations. In the case of blazars inclined at relatively 
large angles such as Cen A, we predict the 7-ray spectra 
extending through the TeV energy range even at large in- 
clination angles provided that the blob velocities are mildly 
relativistic (as observed in the case of Cen A). For slowly 
moving blobs, the GeV 7-ray emission at large inclination 
angles can be even larger than for small angles. However, 
the TeV 7-ray emission behaves in the opposite way. There- 
fore, the model predicts even higher chances of detection of 
nearby blazars with blobs moving with mild velocities and 
inclined at large angles in respect to those once inclined at 
small angles. 

Assuming that the injection rate of primary electrons 
within the blob is constant during the propagation within 
the jet, the 7-ray flares last for a few hours up to a few 
days depending on the range of injection distances from the 
accretion disk. However, the 7-ray spectrum does not change 
significantly with the development of the flare. It is expected 
that at large inclination angles the 7-ray flares observed at 
the GeV energies should precede the 7-ray flares at TeV 
energies. This is caused by the strong absorption of TeV 7- 
rays propagating at large inclination angles at the beginning 



of the flare (when the relativistic electrons in the blob are 
still relatively close to the accretion disk). 

In the case of blazars observed at relatively small in- 
clination angles and with fast moving blobs (e.g 3C 279), 
the 7-ray emission features differ significantly. Such blazars 
are typically located at cosmological distances. So then, 
additional effects are caused by the propagation of 7-rays 
through the Extragalactic Background Light. Let us at first 
concentrate on the 7-ray emission features in the source 
frame (before propagating in the EBL). The 7-ray emission 
during the initial stage has similar intensity (the plateau 
phase) which is followed by an exponential decay whose time 
scale depends on the range of injection distances along the 
jet and the velocity of the blob. The 7-ray spectrum has 
similar shape during the plateau phase extending to larger 
energies with the development of the phase. However, dur- 
ing the exponential decay of the flare, the 7-ray spectrum 
becomes steeper and cut-offs at lower energies. Therefore 
we predict significantly longer 7-ray flares at the GeV en- 
ergies than at the TeV energies. The 7-ray spectra in the 
observer's reference frame (after propagation in the EBL) 
are drastically absorbed at energies above a few hundred 
GeV. Therefore, the shapes of the 7-ray spectra at different 
stages of the flare become more similar to each other. 

In order to envisage whether the considered model pa- 
rameters are reasonable, keeping in mind the limitations of 
the simplicity of the scenario, we compare our results with 
the spectral features of two blazars visible at a small and 
large inclination angles, Cen A and 3C 279. The GeV- TeV 
7-ray spectrum observed from Cen A can be well described 
by our model in the case of the observer located at the angle 
< 40° for the electron injection spectrum with the spectral 
index —3 at the range of distances between 1 — 300 ri n from 
the base of the jet (see Fig. 3). Note that the estimates of 
the inclination angle of the jet in this source lay in the range 
10° — 80°. So then, our calculations are consistent with the 
lower values of the inclination angle. The spectra calculated 
for larger angles are clearly too steep. Unfortunately, Cen 
A is a relatively weak source. Therefore, the 7-ray spec- 
tral information at a flare state and the light curves are at 
present unavailable. We have to wait for the next genera- 
tion Cherenkov instrument (e.g. CTA) and/or a strong flare 
from this object in order to test some predictions of our 
model mentioned above (e.g. hard to soft evolution of the 
high energy 7-ray spectra with developing of the flare). 

We also show satisfactory fit to the GeV- TeV 7-ray 
spectrum reported by the EGRET and the MAGIC tele- 
scopes from 3C 279 during the flare state and predict how 
these spectra should evolve in time (see Fig. 8). The spectral 
information from distant blazars (such as 3C 279) is partially 
hidden by the propagation effects in the EBL. Therefore, it 
will be more difficult to confront predictions for the spectra 
(e.g. its evolution in time) and light curves (e.g. the plateau 
and tail emission) in the 7-ray range expected in terms of 
our model with the observations. 

In the case of Cen A, we modeled the GeV- TeV 7- 
ray spectrum in terms of the cascade model initiated by 
primary electrons (pure external Compton e pair cascade 
model, this work) and the cascade model initiated by pri- 
ma ry 7-rays (produced in the jet, e.g. in the SSC model, 
sot' ISitarek fc Bednarek! (|20ld )). There aren't any definitive 
differences between these two scenarios. Both can describe 
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the spectral featu res equally well (compare F ig. 3 in this pa- 
per with Fig. 4 in lSitarek fe Bednarekl (^Old ) - ). However, the 
injection spectrum of primary particles have to differ signif- 
icantly. Moreover, it is expected that the application of the 
SSC model as a source of primary 7-rays for considered here 
cascade scenario may have additional consequences. For ex- 
ample, SSC scenario may produce additional soft radiation 
field which may compete with the accretion disk radiation. 
As a result, other component in the broad band high en- 
ergy spectrum might appear and some effects on the IC e 
cascade spectrum produced only in the disk radiation is ex- 
pected due to an additional cooling mechanism of secondary 
leptons. T hey have not been included in the cascade calcu- 
lations bv lSitarek fe Bednarekl (|2010h . 

We have discussed the 7-ray emission features in terms 
of the anisotropic, time dependent, IC e pair cascade model 
in the case of a simple, constant lepton injection for a range 
of parameters. The 7-ray light curves can change signifi- 
cantly in the case of a more complicated model which in- 
cludes, e.g. the change of the injection rate of electrons along 
the jet, the extended injection of electrons at the base of 
the jet, or the escape of electrons from the blob. Moreover, 
blobs can propagate along the jet in a more complicated 
manner than considered in this simple scenario. For exam- 
ple, the blob may decelerate during its propagation along 
the jet and its viewing angle may also change as postulated 
in some works. Therefore, detailed modeling of the devel- 
opment of 7-ray flares from specific blazars in such a more 
complicated version of the external Compton model have to 
wait for a more complete knowledge on the processes within 
the jets of AGNs. 

We would like to stress an importance of the simulta- 
neous multiwavelength observations of blazars observable at 
different inclination angles. They would allow a more de- 
tailed modeling of the processes responsible for the pro- 
duction of 7-rays in blazars. Simultaneous GeV and TeV 
observations of flares from blazars can efficiently limit the 
acceleration region and the dynamics of the injection rate. 
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APPENDIX A: IC e± PAIR CASCADE IN THE 
DISK RADIATION 

We consider the process of interaction of electrons injected 
close to the jet base with the soft thermal photons emitted 
from the surface of the accretion disk. Due to the large opti- 
cal depths, electrons interact efficiently with soft photons in 
the inverse Compton process producing 7-rays. These 7-rays 
again interact with the same disk radiation producing next 
generation of leptons. As a result, the IC e pair cascade 
develops in the anisotropic radiation field created by an ac- 
cretion disk in the whole volume above the disk, i.e. outside 
the jet volume. Here we show the method which is used to 
track the development of particles in such a cascade. The 
method allows us to produce 7-ray spectra seen by external 
observer at arbitrary direction with respect to the accretion 
disk. 

Al e pairs from 7 — 7 interaction 

Let's consider a 7-ray photon with energy, E^, propagating 
at an arbitrary direction at a given distance from the ac- 
cretion disk (see Fig |Al|) . Differential energy density of soft 
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Figure Al. Schematic presentation of the geometry of inter- 
action of the 7-ray photon with the soft photon with energy e 
coming from the accretion disk. Soft thermal photon is emitted 
from the unit surface dS c at the angle, 0' , to the disk surface. The 
emission place is defined by the distance, r, from the center of the 
accretion disk and the angle, ip. The location of the 7-ray photon 
is defined by the distance from the axis of the accretion disk R 
and the height above the disk surface H. The distance between 
dS e and the location of 7-ray photon is marked by R! . Photons 
interact at the angle 9. 



photons from the disk, seen by the 7-ray photon, is 



dn 



cos 6' 



'" dSededV h :i c 3 exp(e/fcT) - 1 R' 2 ' ^ 

where h is the Planck constant, k is the Boltzmann constant 
c is the velocity of light, T = T(r) is the temperature profile 
on the disk surface at the distance, r, from its center (see 
Sect. [21), and other parameters are defined in Fig. EH 

At first, we determine the propagation distance, Lint, 
after which, a 7-ray photon with energy _E 7 interacts with 
the soft photon from the disk, 

\{x)' 1 dx = - ln(l - P r ), (A2) 



/ 

Jo 



where P r is the random number, and X(x) is the mean free 
path for e pair production in 7-7 collision at the propa- 
gation distance, x. X(x)^ 1 is calculated following by general 
formula, 



1 (1 — cos9)ndedip rdr, (A3) 



where cr 77 is the cross section fo r 7 — 7 — > e + e pair pro- 
duction (| J audi fc Rohrlichllilfflil ). and rin and r out are the 
inner and outer radius of the accretion disk. \(x)~ 1 can be 
re-written into the form suitable for numerical calculations, 

A " = 8^ [^drr {k Tffjd,I {m - C osO)^ 

(A4) 

where <tt is the Thomson cross section. In the above for- 
mula the internal integration over the energy e was sub- 
stituted by integration over resulting electron velocity /? 
in the center of mass system. The value of internal in- 
tegral 1(b) is parametrized with a single parameter b = 
2m 2 /[E 1 {l — cos#)fcT], m c is the electron rest energy: 



P 



■p 2 Y exp(6/(l - P 2 )) 



-g(P)dp, (A5) 



g(p) = (1 - p 2 ) 2p (p 2 - 2) + (3 - p 4 ) In 



1+1 
1-P 



(A6) 



If Li n t is within the simulation area (defined by the maxi- 
mum distance 10 4 ri n ), the 7-ray is absorbed in the e ± pair 
production process, otherwise it escapes the accretion disk 
radiation field. 

We determine the energies of secondary electron and 
positron in the case of absorption of 7-ray photon. At first, 
we draw the place on the disk from which interacting soft 
photon originates and afterwards its energy. The probability 
of an interaction of 7-ray photon with soft photons coming 
from a given location on the disk is proportional to the func- 
tion under integral in Eq. (|A4|) . 



-drdtp = f~ n (r,ip)drdtp. (A7) 



dP <xrT 3 1(b) (1- cos <9)^!p 

We use the rejection method in order to determine the pa- 
rameters of the soft photon emission place on the disk i.e., 
r and p. To do that, we estimate the maximum value of 
/ 77 . Using the known temperature profile of the disk sur- 
face and numerically determined maximum possible value 
for 1(b) < 2.4, we get the condition / 77 (r, (p) < / max = 
2 • 2.4 • Tf n r 2 n / H 2 . Afterwards we draw random values of r, 
tp and P r until the condition / 77 (r, cp) < P r ■ / ma x is fulfilled. 

In the next step, we determine the velocity, P, of leptons 
produced in the center of mass system. It is directly linked 
to the energy of interacting soft photon emitted from the 
specific place on the disk by the formula e = bkT/(l — p 2 ). 
Distribution of p is proportional to the function under inte- 
gral in Eq. lA5l We determine the maximum value of, /2, m ax, 
for fixed value of parameter b. Now, we randomly sample the 
velocity of lepton produced in the center of mass system. 

In order to determine energies of secondary leptons, at 
first the scattering angle (in the center of mass system), 
#sc (cos# sc = /Use), of the produced lepton has to be found. 
/j BC is randomly sampled from the angular distribution (see 
iJauch fc Rohriichl (| 19761 )1, 

da ^ p(l~P^j c + 2p 2 (l-p 2 )(l-fi 2 c )) ^ 



(1-P 2 A) 2 



by reversal method of the above probability distribution. For 
known P, the energy E' and momentum p' of lepton in the 
center of mass system is determined. The velocity of the cen- 
ter of mass system, /3 C ma, and its Lorentz factor 7 C ms, is cal- 
culated from p cm3 = y/E 2 + e 2 + 2E 1 e cos 9/(E~, + e). This 
allows us to determine energies of leptons in the accretion 
disk reference frame by applying obvious transformation, 



E„ 



r y C ms(E' ± PcmefJ-scP)- 



(A9) 



It is assumed that directions of secondary pairs produced 
in 7 — 7 interaction are the same as the original direction of 
the primary 7-ray photon. 



A2 7-rays from IC scattering 

The method of simulation of the parameters of the secondary 
7-rays produced by leptons in the inverse Compton Scatter- 
ing process is essentially the same as described above for the 
production of e pair in 7 — 7 collision. Let's consider a lep- 
ton with a Lorentz factor, 7 (and the corresponding velocity 
P) interacting with disk thermal photon. The energy of pho- 
ton in the lepton's frame before and after scattering is equal 
to e = u}'m e and e' ac = oj' ac m e respectively, m e = 511 keV. 
bJ is related to the photon energy e in the disk frame by a 



Time dependent gamma-ray production in the IC e^ pair cascade 17 



transformation: u' — e' /m c — 7e(l — f3/j,)/m c = ae, fi is a 
cosine of the angle between direction of the soft photon and 
the lepton. 

We use general Eq. (|A3|) for calculating the mean free 
path, but with the cross section replaced by the to t al cro ss 
section for the IC process (see IJauch fc Rohrlichl (| 19761 )). 
For numerical calculations, we use the following formula for 
the IC mean free path, 

Afcl = 2(Sp / rf ( lnr )r 2 | Ml-cos^^/jc, (A10) 
where 



he 



s(oj') 



exp(t<j'/6) 



-duj' ' 



(All) 



b = akT. The total cross section for IC is aic = 2irros(uj'). 
The interaction place of lepton with the soft photon is deter- 
mined according to Eq. (|A2|) . In the next step, we simulate 
the place (r, <p) from which the scattered soft photon has 
been emitted. Again the rejection method is applied to the 
probability function proportional to the function under in- 
tegrals in Eq. (|A10I) . 



dP oc r 2 (l - cos 6») 



cos 9' 



Tied hi rdip. 



(A12) 



After the interaction place is fixed, we simulate the energy 
of interacting soft photon in the lepton's frame by drawing 
to' from the distribution (applying the rejection method), 



dP oc s(uj') 



exp(o//6) 



-dw' . 



(A13) 



Afterwards the scattering angle of the soft photon in the 
electron reference frame is evaluated from the distribution, 



da M sc + w - /J, sc u + 
oc — 



1+u 



(A14) 



d/isc (1 + Ci/ — /i sc Cj') 2 

where fi BC = cos(9 sc . We rewrite Eq. (|A14jl in a variable 
x = 1 + u' — /i S c<^>' and apply the reversal of the cumulative 
probability distribution method, 

da d[i 3C ^ _ 

~ 2io ,i x 2 



D[x] 



._ = — ^— ( 2(u/ 2 - 2J - 2)x 2 ln(x-) 
+ ((x - l)(4u'x + u/ 2 (l + x) + 2x(l + a;))) . (A15) 



The value of x (and related to it scattering angle fi ac ) is ob- 
tained by solving the following equation for random number 

Pr, 



D[x] = P r ■ D[l + 2w'] 



(A16) 



Now, we are able to calculate the energy of 7-ray photon 
produced in the IC process in the disk frame from, 



Ey = 7 £ L(1 + /Jcosq) 



(A17) 



where the energy of the scattered soft photon in the electron 
rest frame is e' gc = e'/[l + e'(l — /i sc )/m c ]. The angle between 
directions of relativistic boost and scattered photon is cal- 
culated as cos a = cos 9 SC cos 9' + sin 9 ac sin 6' cos ip a c (see 
Fig. IA2[I . The interaction angle in the electron rest frame is 
related to the corresponding angle in the disk frame 9 with 
the transformation: cos9' — cp' x /e' = 7e(cos6 ) — /3)/e'. 




Figure A2. IC scattering as seen in the electron rest frame. The 
direction of the relativistic boost (the electron movement direc- 
tion in the disk frame) is shown with K e , 7-ray direction is K!y. 
The disk photon is interacting at the angle 0' to the boost direc- 
tion. After scattering its direction is determined by the angles 9 SC 
and ip S c- The scattered photon is inclined at the angle a to the 
boost direction. 



APPENDIX B: TIME DEVELOPMENT OF THE 
IC E ± PAIR CASCADE 

In order to obtain the gamma-ray light curves and the time 
dependent spectra we have to track the time development 
of the cascade above the accretion disk. Each particle has 
an assigned clock, t, running in the reference frame of the 
accretion disk. The injection process occurs continuously in 
a blob moving with the velocity v along the jet in the range 
of distances between Hi and Hi. The injection time of elec- 
trons at the height, Ho, is equal to 



t = (Ho - HJ/v. 



(Bl) 



If a new particle is created in the cascade, its clock is set to 
the time of its parent particle during the interaction. 

In order to switch from particle clock time, to a time of 
external observer we have to apply a correction for the emis- 
sion distance and direction. In reality the source is observed 
under a specified inclination angle £o- In order to obtain 
enough statistics, the 7-ray spectra and the light curves are 
obtained by summation over a small range of cosine of ob- 
servation angles (e.g. A(cos£) = 0.1). Within such a small 
A(cos^) interval, the time and spectral characteristics of 
the 7-ray emission does not change significantly. Below we 
describe the method which has been used for the correct ar- 
rival time evaluation in order not to introduce any artificial 
delays due to integration over a range of £ angles. 

In case of the lack of the magnetic field in the region 
above the disk, the secondary leptons propagate exactly as 
7-rays i.e., along straight lines. Therefore, each cascade ini- 
tiated by a primary electron in the jet propagates along a 
straight line starting from a point on a jet. The extrapolated 
height above the accretion disk H e can be calculated from 
(see Fig. lBlj) : 



H e = H — dcos£, 
d = i? sin -1 £. 



(B2) 
(B3) 



The gamma-ray A, produced at the time t, is equivalent to 
a gamma-ray sent from the point P (at the height H e above 
the disk) at the time t' = t — d/c. 

After extrapolation of directions of escaping 7-rays to 
the locations in the jet, we calculate a second correction for 
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Figure Bl. Position dependent time correction applied in order 
to obtain the 7-ray light curves at the observer reference frame. 
Two 7-ray photons, A and B (thick arrows), are produced in 
two sub-cascades at similar inclination angles and at a distance, 
d, from the location of the parent electron (H e above the disk 
surface). The 7-ray marked by, A, is characterized by the height 
over the accretion disk H, the distance from the jet axis R, and 
the inclination angle £. The inclination angle £0 describes the 
center of the cosine angle interval in which 7-rays escaping to the 
observer are collected. 

the height of H e . Instead of using the inclination angle of 
the individual 7-rays, we use the central value, £0, of a range 
of the inclination angles A(cos£). In this way, both photons 
A and B will be corrected by the same time difference of 
.ff e cos£o- Therefore the final correction which is applied to 
all 7-ray photons is 




(B4) 



